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This investigation deals with some exact solutions of the equations governing the steady plane 
motions of an incompressible third grade fluid by using complex variables and complex functions. 
Some of the solutions admit, as particular cases, all the solutions of Moro et al. 
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I. INTRODUCTION 

The governing equations that describe flows of New- 
tonian fluids are the Navier- Stokes equations. The me- 
chanical behavior of many real fluids, especially those of 
low molecular weight, appears to be accurately described 
by these equations over a wide range of circumstances. 
There are, however, many real substances which are ca- 
pable of flowing but which are not at all well described 
by the Navier- Stokes theory. Due to this reason many 
fluid models have been proposed and studied by differ- 
ent authors. Among these, the fluids of differential type 
have received much attention. Interesting studies of dif- 
ferential type fluids are given by Rajagopal Erdogan 
@], Rajagopal and Gupta @], Bandelli Siddiqui and 
Kaloni [5[, Benharbit and Siddiqui @, Ariel B[8L Fete- 
cau and Fetecau [9] and Hayat et al. [ToL [TlT fl2|]. The 
fluids of third grade, which form a subclass of the flu- 
ids of differential type, have been successfully studied in 
various types of motions. 

Moro et al. [l3[ determined some exact solutions of 
the equations governing the steady plane motion of an in- 
compressible third grade fluid employing hodograph and 
Legendre transformations. 

The aim of this paper is to present some exact solu- 
tions for steady plane motions of incompressible third 
grade fluids. For this, the complex variables and com- 
plex functions are used. Some previous solutions can be 
obtained as special cases of our solutions. In sections 2 
and 3 are presented the constitutive equations and the 
expression of the vorticity function uj. Section 4 contains 
solutions and their graphical illustrations. 



II. CONSTITUTIVE AND GOVERNING 
EQUATIONS 

The fluids of grade n, introduced by Rivlin and Er- 
icksen LL4L are the fluids for which the stress tensor is 



a polynomial of degree n in the first n Rivlin-Ericksen 
tensors defined recursively by 



A n 



(diUj + djUi)ij ; 
d 



dt 



A n - 



A n -iL + Zr A n _i , n > 1 (1) 



where ^ = dt + u ■ V denotes the material derivative and 

L = {djUi)i^j ; if = (diUj)i j j . 

Physical considerations were taken into account by Fos- 
dick and Rajagopal [15[ in order to obtain the following 
form of constitutive equation for the Cauchy stress T in 
a third grade fluid: 



T = -pi + Mi + ai A 2 + a 2 A\ + /M 3 • 
+HAiA 2 + A 2 Ai] + fo(trA\)A 1 



(2) 



where —pi is the spherical stress due to the constraint of 
incompressibility, \i is the dynamic viscosity, ai(i = 1,2) 
and — 1,2,3) are the material constants and the 
three Rivlin-Ericksen tensors Ai, A 2 and As are given in 
equation (1). 

Furthermore, a complete thermodynamic analysis of the 
constitutive equation (2) has been given by Fosdick and 
Rajagopal [15[. The Clausius-Duhem inequality and the 
assumption that the Helmholtz free energy is a minimum 
in equilibrium provide the following restrictions, 



(3) 



/i>0, ai>0, \ai +a 2 | < V 24 ^, 
ft = 02 = 0, 03 >0. 

Thus equation (2) becomes 

T = —pi + ol x A 2 + a 2 A\ + [fj, + /3 3 (trAl)]Ai 



where ji e ff = \i + /3s(trA\) is the effective shear- 
dependent viscosity. 

The velocity field corresponding to the motion is given 
as: 
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v = v(x, y) = 00, y), v(x, y), 0) . 
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III. FLOW EQUATIONS 

The basic equations governing the steady plane motion 
of a homogeneous incompressible fluid of third grade, in 
the absence of body forces are [l3j] 
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where 
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(5) 



(6) 



M = 8 
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Let 



z = x + ty and 
Then from Stallybrass 
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Equation (9), on utilizing equation (10), becomes 



where 
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^-a 1 (u\7 2 u + v\7 2 v)-(3a 1 + 2a 2 )^-+p (7) 



dzdz 
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M = 32 



<9 2 V> <9 2 ^ 
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Equations (5) and (6) are non-linear partial differential 
equations for three unknowns u, v and p as functions of 
x and In equations (5), (6) and (7), the viscosity p 
and the material constants ai, 0^2, A : A : A satisfy the 
constraints given in equation (3) [15L fl6|. 
Equation (4) implies the existence of a stream function 
ip(x,y) such that 



dip dip 
— and v = - — 
oy ox 



(8) 



Equations (4) and (5), on utilizing equation (8) and the 
compatibility condition ~ £ = i 1 , yield 



f dip du dip du \ 
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-0 (9) 
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IV. SOLUTIONS 

On integrating equation (12), we obtain the general 
form of the exact solutions 



' dz dz + A + A 



(13) 



where A and A are complex functions and A is the com- 
plex conjugate of A. 

To determine the solution of equation (11), our 
strategy will be to specify the vorticity function u and 
determine the condition, which the functions A and A 
must satisfy. The condition is obtained from equation 
(11) utilizing equation (13). 

I. When vorticity u is constant say u Q , then equation 
(13) yields 



ip — —-u zz + A + A 
Equation (11), on utilizing equation (14), yields 



(14) 



p3 



/d*A d 4 A\ 
\~dz^ + l?z T ) 



dz 2 



dz 2 
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When p 3 ^ 0, then 

lciiz 4 a 2 z a a 3 z^ . , . 

A= ^ + ^ + ^-+ a 4^ + a 5 (15) 

where a^s are all complex constants. The velocity com- 
ponents u and v are given by 



ia\z 3 a 2 z 2 



lcliz 3 a 2 z 2 



6 3 

a 3 z - a 4 f ; 



• a 3 z + a 4 



uo Q , _ x rmi^ 3 a 2 £ 2 
v = —{z + z) - I — h — ^— + a 3 2 + a 4 



miz 3 a 2 £ 2 



6 3 

ajz + «4 1 



6 3 

Equation (14), on utilizing equation (15), yields 

ujozz i(aiz 4 — a\z A ) (a 2 z 3 + a 2 z 3 ) 



(a 3 z 2 + a 3 z 2 ) 



24 



+ (a4^ + a 4 z) + a 
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(16) 



II. When uo is non-constant, the solutions of equation 
(11) are determined as follows: 

(i) When uo = m\z + ralz , the equation (11) 
yields 

_ r m?z 2 cM^i r . . 

im j h mi j> = A 1 6m i mi (mi z + m^j — 

< 

mi - 



/ a 3 A#i <9 3 A<9 2 A\ 



+64 



d 3 A d 3 A 

si 3 " si 3 " 



(miz + miz)| 



where A 



(33 



and mi is the complex constant. 



The L.H.S of equation (17) suggests to assume 
dA 



Xiz 2 + X 2 z + A 3 



<9z 



This on putting in equation (17), gives 



(17) 



where a = a*, + 05 . 

The stream function -0 in equation (16) is represented 
graphically in Fig.l with uoq = —1, a\ = 1 + 2*,, a 2 = 
1 + ^, «3 = 1 + 5*,, a 4 = 2 + 0.5*,, a = 2 and x , 2/ G [—1, 1]. 



Ai 




Figure No. 1: Graphical representation of equation (16) 

When /3s = 0, then equation (15) is identically satisfied 
and the complex A becomes arbitrary, due to which A 
enables us to construct a large number of streamfunction 
ip and hence a large number of solutions to the flow equa- 
tions. We mention that by taking uo = or appropriately 
choosing complex constants a 3 , as, and a\ = = a 2 
in equation (15), or by taking A = talogz (for /? 3 = 0), 
we get all the solutions of Moro et al. 131] . 

Further more, if we take the function A = (ci + tc 2 )z 2 
or A = (ci + *,c 2 )lnz and choose appropriately the 
constants or apply the appropriate boundary conditions 
we get the plane Couette flow, the flow due to a spiral 
vortex at the origin and the flows having streamlines 
as a family of ellipses, concentric circles, rectangular 
hyperbolae. 



8mi 



A 2 = 40*,Ara?, A 3 = 0. 



The solution of equation (17), therefore, is 

9 3 

A = - 20i\m 2 z 2 + m 2 

24mi 

where m 2 is a complex constant. 
Equation (13), using (18) becomes 

\m\z 3 + m\z 3 \ 

i\) = - — (rmz + rnTz) — — 

8 24mimi 



+2()L\(m 2 z 2 - mlz 2 ^ 



(18) 



(19) 



where m = m 2 + m^ . 

The stream function ip in equation (19) is represented 
graphically in Fig. 2 with mi = 1 + 2l , m=l, A = 0.3 
and x : y G [—1, 1]. 




-1 T -1 

Figure No. 2: Graphical representation of equation (19). 
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(ii) For uj = B(z + z) , the constant B being real, the 
equation (11) becomes 



Im 



c Bz 2 dA-\ 
[- — + — > = A 



-8B( 
+64 ( 



8 02J 
d 2 A d 2 A 



{6B 2 (z + z) 



dz 2 dz 2 
0+40^4 d 2 A d 3 A 

dz 2 dz 3 dz 2 dz 3 . 



\ / d 3 A d 3 A\ 



. d 3 A d 3 A 



2 d 4 A , — d 4 A\ 



+32 ( 



dz 3 dz 3 
d 2 Ad 4 A d 2 A d 4 A 



&z 2 dz 4 dz 2 0(2) 
64 /0 2 y4 2 4 A d 2 A 2 d 4 A 
~B\ dz 2 dz 4 + dz 2 Ih 4 



z) 4 ) 

)} 



(20) 



The L.H.S of equation (20) suggests ^ to be a polyno- 
mial in ~z of degree two and therefore on substituting 



dA 



hz 2 + l 2 z + A 4 , 



in equation (20), we get 
B 



h = , 

8 ' 



Hence 



— Bz 3 

A = — — + 20l\B 2 z + ra 3 



A 4 = 



(21) 



where ms is an arbitrary complex constant. 
Using equation (21), equation (13) implies, 



Bzz. 



(z+z) - — (V + z 3 J + 20iB 2 \(z-z)+n (22) 
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where n = ms + ms . 

The stream function ip in equation (22) is represented 
graphically in Fig. 3 with B — — 2 , A = 2, n = 1 and 

x e [-10, io] , ye [o, io] . 




Figure No. 3: Graphical representation of equation (22). 
( hi ) When u = D(z + z + E), DandE being real, 



then equation (11) yields 

Dz 2 DEz dA 



Im\ 



-8D 



'd 2 A d 2 A\ 16D ( z d 3 A _d 3 A\ 

- dz 2 dz 2 ) V dz 3 dz 3 



0+40+4 0+40+4\ / 0+40+4 
dz 2 dz 3 + dz 2 dz 3 ) +Qi { z dz 3 dz 3 



+64 ( 

d 3 Ad 3 A\ ,„ „ „.„0 3 y4 3 A 



+2 



dz 3 dz' 



+ 4L> 2 £' + 64.fr 

>2 3 / 



0^ 3 02 3 



A ^ ( _d 4 A 2 d 4 A\ QA( d 2 A 2 d 4 A 

V_ ~&¥ + Z dz 4 )_ D \~dz^ffz T + 

2 A 2 d 4 A\ ( d 2 A 2 3 4 A < 2 A 2 4 A ^ 
+ ^_<> ^ 7 J + 1 



02 2 0z+ 



02 2 0z 4 0z 2 0Z 4 



Following the same procedure as that of previous case we 
find 

_ nT3 DEz? 

A = — — + 20l\D 2 z 2 + 20aL> 2 £z + ra 4 (23) 

24 8 

where m 4 is an arbitrary complex constant. 
Using (23), equation (13) becomes 



+. 



Dzz 



(z + 2+2£)-g(z 3 + z3) 



DE 



-f z 2 ) + 20l\D 2 (z 2 - z 2 ) + 20l\D 2 E(z -z) + q (24) 
where g = m 4 + raj . 

The stream function ip in equation (24) is represented 
graphically in Fig. 4 with D = E = 1 . A = 2 , q = — 1 
and x, 2/ G [-10, 10] . 




Figure No. 4: Graphical representation of equation (24). 

(iv) When ui = Bl(z — z), then equation (11) is satisfied 
provided 

A = -^-20i\B 2 z^ + m 5 

where m$ is an arbitrary complex constant and the 
streamfunction ip is given by 



Bt 



(z-z)+^(z 3 -^)+20l\B 2 ( 



-zz 



z 3 -z 3 )+r 



(25) 
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where r = m 5 + . 

The stream function ip in equation (25) is represented 
graphically in Fig. 5 with B = — 5 , A = 3, r = 10 and 
x G [-2, 10] , 1/6 [0, 8] . 




Figure No. 5: Graphical representation of equation (25). 

( v ) For uo = Bln(zz)+Di , the equation (11) is satisfied 
provided 

A = uiqItiz + rrij (26) 

where rriQ, B, D±, are real constants and mj is complex 
arbitrary constant. 

According to equation (26), the streamfunction ip be- 
comes 

ip = — ^zz^ln(zz) — 2 j — —^-zz + rriQln(zz) + s (27) 
where s = rrij + m7 . 

The stream function -0 in equation (27) is represented 
graphically in Fig. 6 with ttiq = 2 , s = 4 and x G 
[-10, 2], 2/G [-1, 2]. 




-0.5 



Figure No. 6: Graphical representation of equation (27). 

The stream function ip in equation (27) represents the 
motion whose streamlines are concentric circles with 



constant speed along each streamline for the appropriate 
choice of the constants. 

( vi ) When uj = Bzz, B being real constant, then 
equation (11) is satisfied provided 

A=—lnz + m 8 (28) 
9 

where rug is a complex constant. 

According to equation (28), the stream function ip be- 
comes 

i> = -^) 2 --ln(=)+t (29) 
lb p \zJ 

where t = nig + Wig . 

The stream function ip in equation (29) is represented 
graphically in Fig. 7 with B = 1, /i=12, p = 1 , t = 2 
and x , y G [—1, 10] . 




Figure No. 7: Graphical representation of equation (29). 



V. CONCLUSIONS 

In this paper, we reconsidered the flow equations of 
Moro et al. [13| with the objective of determining some 
exact solutions. 

For this purpose the vorticity function uj and the 
stream function ip are expressed in terms of complex 
variables and complex function. The condition which the 
complex functions must satisfy is determined through the 
equations for the generalized energy function h by using 
the compatibility condition q x q v = qJq x • 

Some exact solutions to the flow equations are deter- 
mined using the condition for the complex functions. The 
solutions presented in this paper admit, as particular 
cases, all the solutions of Moro et al. [l3[ by appro- 
priately choosing the complex functions or the arbitrary 
constants therein. 
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